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ABSTRACT 
A two-dimensional dynamic model of a towed buoy designed to remain 
submerged just below the water surface is developed. The study includes 
the model development, the analysis of certain model parameters, the 


development of a tension force, and the range of stability of the model to 


initial condition perturbations. 
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mass moment of inertia about Y axis, in slug: ft? 


angular acceleration in rad/sec* 
drag, in pounds 


coefficient of lift 

coefficient of drag 

coefficient of drag of the antenna 

coefficient of the pitching moment about the towpoint 


mass density of water, approximately equal to 1.99 
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mean chord length of the buoy 
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I, INTRODUCTION 


The purpose of this investigation is to attempt to mathematically 
model and computer-simulate the two-dimensional motion of a submarine- 
towed buoy. This has a very specific military application in submarine 
communications systems; the basic system configuration aS a submarine 
communication link is illustrated in Fig. 1. It is desired that the buoy 
stay underwater at all times, and remain within 1 ft. of four feet below 
the water surface; this gives the buoy an operating range from three to 
five feet below the surface. An additional system requirement is that the 
buoy system stay within its allowable operating zone while following 
waves up to sea-state four. 

Many investigators, such as Lacey, have dealt with ship-towed or 
airplane-towed devices, where the buoy is the "bottom-following" type. 
Pode and other early investigators studied the steady-state configuration 
of cables in a steady uniform stream; and more recently Schram developed, 
by the method of characteristics, a two-dimensional dynamic cable-buoy 
model. However, he treated the buoy as a vector force and did not go 
into the actual buoy dynamics. 

Steady-state studies do not show the response of the buoy and cable 
to perturbations such as sea state. It is this response that will prove 
most useful in the study of buoy and cable parameters. This study has 
been primarily concerned with developing the dynamic mode! of the buoy; 


and the tension force generated by the cable has been loosely approximated. 
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The model is of a buoy presently being developed for the U. S. Navy; 
thus, the dimensions and other data were obtained from an existing 
prototype. 

This study develops the two-dimensional model around two basic 
equations; force equals mass times acceleration, and the moment equals 
the mass moment of inertia times the angular acceleration. The model 
block diagram was built around these two equations, several basic 
hydrodynamic relationships, and the wind-tunnel and dimensional data 
supplied on the buoy. Once the block diagram was formulated, it was 
then ready to be computerized. The model was simulated using the IBM 
360-67 at the Naval Postgraduate School, and DSL (Digital Simulation 
Language); DSL is a very simple, but versatile simulation language 
developed by IBM. Once the model was computerized, the lengthy and 
time-consuming task of optimizing certain parameters and finding a real- 
mystic formulation of the tension force followed. The determination’or the 
tension force and its coupling with the buoy has been one of the most 
difficult problems in the formulation of the buoy-cable model. This has 
been a problem in all buoy-cable models where the dynamic response is 
of interest; this is mainly due to the fact that it is not exactly known how 
the forces generated by the buoy and cable react on and with each other. 
The development of the model and the results of the simulation are 


presented. 


13 


II, EQUATIONS OF MOTION AND PROBLEM SETUP 
The purpose of the buoy is to remain four feet below the water 

surface so that a protruding antenna may maintain two-way communica- 
tions; when the buoy iS operating in high seas it must follow the sea 
surface. The body, as Shown in Fig. 2., is basically a fat wing that 
produces dynamic lift according to its angle of attack. This angle of 
attack may be changed by torque produced by the fin. Besides the 
dynamic lift, produced by the body angle of attack, there is the static 
lift of positive buoyancy. The buoy model is based on the two basic 
equations: 

F=ma 

M=l, Q 
The first equation is used to describe movement in the horizontal and 
vertical planes, the X and Z planes respectively. The second is used 
tordescribesthe rotation of the»body about the Y-axis; these axis-are 
shown in Fig..2s The three hydrodynamic forces of lift, drag, and moment 
are defined as follows: 


= 2 
L= $C, pv’S 
. 4 2 
D=$CppVv's 
M= $C, pv?sC 


Along with the tension force produced by the cable and the buoyancy of 


the body, these three equations describe all the forces acting on the 
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body. Before going into the complete buoy block diagram, the body 
coefficients will be put into a more convenient form, and the fin servo 


loop, which controls the buoy attitude, will be examined. 


A. BODY @@EFFICIENTS 

In both aerodynamics and hydrodynamics, there are various dynamic 
coefficients; these are determined by the shape of the body involved. 
There are tables of coefficients for certain common wing shapes. However, 
when an odd body configuration is being studied, it may not be possible 
to apply the standard shapes in order to derive overall coefficients for 
the body. In this case the coefficients of the body must be determined 
by wind-tunnel tests. 

This was the case with the system being studied; the body is of 
non-standard shape, but the fin is the standard NACA 0012 shape. The 
waree cCoeiticients, Cy. Ch: and VE were measured at various body 
and fin angles. The coefficient of lift and the coefficient of moment 
include the effect of the positive buoyancy force of 490 pounds. The 
effect of the antenna drag is not included in any of the coefficients. 

The values of these coefficients are shown in Tables I-III as a function 
of body angle of attack and fin angle of attack. It was decided that for 
computer implementation, these coefficients would be easier to handle 
if they were in equation form. This was done by applying a least-squares 
fit twice to each data set; the resulting equations may be found in Sub- 


routine Tink of the digital program of the model. The resultant curves and 
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data points for constant body angle of attack and varying fin angle of 


attack are shown in Figs. 3-5. 


By fw SERV@ "LOOP 

The diagram of the fin and the servo-control system is shown in 
Fig. 6. The purpose of the fin is to orient the body So that it will stay 
four feet below the water surface. The fine produces a moment on the 
buoy which causes the body angle of attack to increase or decrease, 
This increase or decrease in body angle of attack causes the body lift 
to increase or decrease and the body rises or descends as is necessary. 

The fin actuator loop is driven by a Size 23, Kearfott motor. The 
loop is designed to have a rise time of 0.04 sec.; this actuator loop is 
shown in Fig. 7. The servo loop used in the buoy model is considered 
to be a K/s* plant with tachometer and unity feedback; this is shown in 


Fig. 8. The system transfer function: 


2 
te: ee 
i+ GH le K/s*(K,s + 1) 

a 
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s* + K KS + Ky 


and the characteristic equation: 
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+K Ks + = 
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Assuming K = 1000 


Ww Kee = 10.00 
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THE K/S* SERVO SYSTEM 


Figure 8 
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Taking a value of € = 0.7: 


2(0.7)(31.62)/1000 


ot 


.0443 
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With Ky and K, set at the above values, the rise time of the servo 


loop is about 0.1 sec. This is not as fast as desired, but or changing 
the two gains the desired rise time may be obtained. Two different 
responses of the fin servo loop are shown in Fig. 9; both have an input 
step of 0.174533 rad. (10 degrees). The first set of gains is used in the 
model simulation. This is because it is doubtful that the inputs into the 
servo loop will require that the fin move that fast; if so, then the buoy 
has probably reached the limits of its surface-following capability. 

The inputs into the fin servo loop are shown in Fig. 6. The fin 
servo loop that was simulated above is composed of the summer, servo 
amplifier, motor and gears, linkage to the fin, and the unity and tacho- 
meter feedback from the fin to the summer. The other three inputs will 
come from gages and electronic sensors in the body. These inputs are 
the depth error, the depth rate, and the body pitch rate. In the model 
Simulation, there are three constants (ADE, ADR, and APR) which multiply 
these three inputs. The values of these three constants are discussed 


later in the study. 
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Figure 9 
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Ill, THEaBwOyY MODEL 

In this development of the buoy model, the effect of water currents, 
the forces generated by the heaving water surface, and the free-surface 
effect are not taken into account. Currents can be considered neglible 
compared to the heaving forces the buoy might experience. The free- 
surface effect is experienced by the buoy because it is a body moving at 
a depth comparable to its length; a more detailed and complicated model 
would have to include this force. The force generated by the heaving 
water surface would not affect the buoy as much as might be expected. 
ims isesecause the body, as gma as its vemical freedonmisme: restricted, 
will act like a particle in the water and will move up and down with the 
Sumace we nus, thesieaving forees will aetually ellp nrarmtam the buoy 
near the water surface, whichis what is desired. If the buoy were being 
towed directly above its towpoint, it would have little or no vertical 
freedom; however, if the buoy is well behind the towpoint, then the cable 
will provide little restriction to the buoy's vertical freedom. Tests on the 
actual buoy have shown that with 475 ft. of cable and the two point ata 
depth of 240 ft., the buoy can heave upto 20 ft. with the waves Without 
the cable decreasing the buoy's vertical freedom of motion by more than 


10 percent.? 


* Soulant, H., Sea Surface Followage by Towed Buoy (Preliminary 
Study), p. 11, March 1969. 
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The buoy model will now be developed; the buoy diagram is again 
produced with some specific dimensions (Fig. 10). Besides the diagrammed 
dimensions, the buoy has the following additional characteristics: 


ne 25. tt° 


it 


S 


C= ono 
Ie = 92,3 7aRit 
mass = 70 slugs 
mass moment of inertia = 975 slug SS 
volume = 17.5 ft® 

The block diagram of the dynamic model is shown in Fig. 11. As 
was previously stated, the buoy model is based on the summation of 
forces in the X and Z planes, and about the Y axis. The summed forces 
in the X and Z directions are divided by the mass of the buoy, and the 
moments about the Y axis are summed and divided by the mass moment 
of inertia. This yields the acceleration in the horizontal and vertical 
directions, and the angular acceleration about the Y axis. The accelera- 
tions are each integrated to give the velocities, and these, when inte- 
grated, give the horizontal and vertical position and the angular 
Cnlentatromy. 

In this first model, the tension force is generated as a function of 
horizontal velocity and body angle of attack. The tension function was 
formed in the Same manner as the coefficient equations, and was based 
on predicted values of the tension in the cable. This tension function may 


also be found in Subroutine Tink of the computer program. 
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The velocities and positions are used to generate the forces and 
moments that are summed in the basic equations. The forces and moments 
on a body are relative to the velocity vector of the medium, and not to the 
— horizontal. Thus, the sum of the squares of the horizontal and 
vertical velocities yields the relative velocity squared, V2. The arc 
tangent of the vertical velocity over the horizontal velocity is taken to 
give the angle of the velocity vector, VAN, with respect to the true 
horizontal. The angular orientation of the body (angle of attack), Al, is 
summed with VAN to produce the angle of attack of the body with respect 
to the velocity vector, AR; ARis in radians. The fin angle of attack, FA, 
is also summed with VAN to yield the fin angle of attack relative to the 
velocity vector, BR. These two angles of attack, AR and BR, are 
multiplied by a conversion factor of 57.29578 deg/rad to yield the angles 
of attack, A and B, indegrees. 

The two angles, A and B, are the arguments which determine the 
coefficient values of body lift, drag, and moment. The coefficients of 
lift and drag are both multiplied times V2 and a constant to produce the 
body lift, BL, and the body drag, BD. The constant is $pS, and is 
egmal to 14.17675. 

ae inement Coeiticient is multiplied times VZ and a consid, *oSC; 
this constant is equal to 84.0725. This body moment, YCM, is Summed 
with the moments produced by the antenna, to yield the total moment, YM, 


acting on the body about the towpoint. 


al 





There is no definite information on the antenna, but itis assumed 
that it will only contribute a drag force; this drag force will contribute to 
the total drag and moment on the body. The antenna coefficient of drag is 
assumed to remain constant at a value of 0.50. Thus, the antenna drag 
is dependent on a constant times the frontal surface area times the 
relative velocity squared. Assuming a frontal width of 2 inches, the 
drag force of the antenna will then vary with the length of antenna under 
water, and V2. The length of antenna will be approximated as the depth 
of the body below the surface. The constant in this case is a(1/6)pC.,, 
which equals 0.04146. The antenna drag is multiplied by the sine and 
cosine of VAN, SV, and CV, to give two moment components. These times 
their respective moment arms yield YADS and YADC. 

The body lift and drag, and antenna drag are multiplied by SV and CV, 
and the resulting six values are summed into the horizontal and vertical 
forces. In this first realization of the model, these horizontal and 
vertical forces are used before being summed with the horizontal and 
vertical components of the tension force, to determine the angle FE with 
which the tension force acts on the body. This means that the tension 
force acts exactly opposite to the vector of body forces; this simplifying 
assumption is used in the initial stages of the model development. 

The depth error and depth rate inputs into the fin servo loop are taken 
from the vertical position and velocity, Zand Zl, respectively. These 
are multiplied by a constant, 0.443, which converts ft and ft/sec into 


psi/ft and psi sec/ft. These are multiplied by gain constants ADE and 


ADR before being summed with the body pitch rate, Yl, times its gain 
constant, APR. These three inputs are fed into the fin servo loop. The 
orientation of the fin causes a torque on the body which either increases 
or decreases the body angle of attack; this enables the body to rise or 
descend. 

This is the initial dynamic model; it changes as more insight and 
knowledge are gained on the model. The first study with the model 


deals with the three coefficients ADE, ADR, and APR. 
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Once the model was digitally programmed and the initial programming 
mistakes ironed out, the first step was to study possible values of the fin 
servo input coefficients. The first computer program realization is shown 
in the program titled "INITIAL BUOY MODEL". The buoy is initiated ata 
depth, ZI, of 4.0 ft, a speed, XII of 20.2536 ft/sec (12 knots), body 
angle of attack, AlI, of 0.10472 rad (6 degrees), and a fin angle of attack, 
FAI, of 0.02 rad (1.15 degrees). The tension force and the coefficients 
of lift, drag, and moment are called out of Subroutine Tink, which is in 


double-precision format. 


A. SERVO INPUT COEFFICIENTS STUDY 

The first attempt in this study was to see if the buoy would stabilize 
with only one input into the fin servo loop; either the depth error, the 
depth rate, or the body pitch rate. This proved fruitless, so different 
combinations of ADE and ADR were tried, and it was found that these two, 
with APR set equal to zero, would stabilize the system. A full range of 
values of ADE and ADR were tried to see which combinations would bring 
the buoy to a steady-state position. A mapping of these values is shown 
in Table IV; the value given for each combination is the depth at which 
the buoy steadied out. The responses for several of these combinations 
are shown in Figs. 12-16. It appears from the depth values that the 
tension force is too strong, Since in every case the buoy steadies out 


deeper than the desired depth of 4.0 ft. This became obvious when the 
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4.3 





4.4 


RESPONSE OF ADE=2.5 AND ADR=0.3 


Figure 16 


steady-state body angle of attack and fin angle of attack were checked; 
both angles were large and showed that the model was working properly, 
but that the buoy was fighting against an unrealistic tension force 
formulation. 

It was decided from the mapping of the two coefficients, ADE and 
ADR, that three sets would be used in further studies of the buoy model. 
The first pair was chosen on the basis that it was on the intersection of 
the values of ADE and ADR which had the widest range of stable values, 
and that this set was towards the center of the mapping. This was 
ADE-1.2 and ADR=0.30. It was then decided that the other two pairs 
would vary in one coefficient away from this point. These other two 
pairs are ADE=1.2 and ADR=0.07, and ADE=2.5 and ADR=0.30. The 
responses with these three pairs of coefficients are shown in Figs. 
14-16. 

The coefficient APR remained to be studied. When this coefficient 
was given a value larger than zero there was no appreciable change in 
the response of the buoy. To check the range of possible values of APR, 
the tension force, T, was changed so that it was equal to the square 
root of the sum of the squares of F! and F2. Thus, the buoy would stay 
at fowr feet below the surface and the effect of APR could be seen im The 
response of the buoy's angle of attack. Both negative and positive values 
of APR were tried, and negative values immediately drove the buoy model 
unstable. Positive values of APR, up to 7.0, were tried with the three 


pairs of depth error and depth rate coefficients. At thiS upper value of 


7.0, the body started to oscillate away from stability about the Y axis; 
it was apparent that this high value of gain was approaching the stability 
limit of the buoy model. Based on this, a median value of 3.0 was 


chosen for APR. The three final sets of coefficients are as follows: 


ADE ADR APR 
dg '2 0 abd 3.0 
hee 2 0.30 3.0 
ye 0.30 3.0 


B, TENSION FORCE FORMULATION 

The next step in this model study was to formulate a new, and hope- 
fully more realistic, tension force, rather than using the predicted 
tension data supplied with the buoy. It was decided that this tension 
function would be a constant plus or minus the distance away from the 
equilibrium point times what will be called a spring constant, SC. The 
basis for this formulation is Shown in Fig. 17. First itis assumed that 
the last portion of the cable connected to the buoy, LO, is straight, and 
that of the total length of the cable this section is the most dynamic. By 
com puting the increase or decrease in ZO and XO, a new cable length, 
LN, may be calculated. If LN is larger than LO then the tension increases 
if it is smaller then the tension decreases. The difference in length 
(LN-LO) is L and this is multiplied times the spring constant, SC, which 
is studied for several different values later on. The new tension force 
now takes the form: 


Parga + Loe 
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THE NEW TENSION FORCE PORMULATION 


Figiipe 7 


The value 1940 is the steady-state tension value that the model 
setiled to in the earlier parameter studies; the angle FSS is the steady- 
state value of FE. This new tension force was implemented into the buoy 
model. The buoy will now be initiated at a depth of zero feet; this is 
possible because no heaving force or surface effect is taken into account 
in the model; other minor changes are made to take this new initial 
condition into account. The new computerized model program is titled 
"STUDY OF T = 1940 + L*¥S5C", The body angle is now being initiated at 
4.8 degrees and the initial fin angle is 0.0 degrees; these two values 
are also the steady-state values for these angles found earlier. The next 
step in the formulation of this tension force was to find an appropriate 
value for the spring constant. The program was run with the three sets 
of coefficients and the new initial conditions, while SC was varied from 
2000. to 0.0 in steps of 250. The two sets of coefficients with ADE=1.2 
were stable for all values down to SC=250, while the set with ADE=2.5 
was Stable down to SC=1750. This third set (ADE=2.5, ADR=0.30, and 
APR=3.0) was dropped from further investigation since this was believed 


to be an unrealistic value for the Spring constant. 


C. STUDY OF INITAL CONDITION PERTURBATIONS 

This final form of the model with the new tension force formulation 
was next examined to see how stable the model was under various initial 
Condition perturbations. It was decided to perturb the initial vertical 
velocity, the body angle of attack, and the fin angle of attack. Three 


different spriny constants were used with the two remaininy sets of 


16 


coefficients; these were SC=500, SC=1000, and SC=1500. It was also 
decided that two different values for the approximation of the dynamic 
length of the cable would be used; these were ZO=100 and ZO=50. (ZO 
Simply because it was used to initiate the model rather than LO). The 
final results were tabulated, and are shown in Table V; the initial angles 
of attack (AlI and FAI) are in degrees, and the initial vertical velocity 
(Z1I) is in ft/sec. 

The responses of the buoy model to some of these initial condition 
perturbations are shown in Appendix A. The twenty-four plots include 
all the maximum perturbations listed for SC=500 and SC=1500. These 
plots do not really tell much about the buoy model. They do show that 
the buoy will settle to steady state fairly rapidly, and that the buoy will 
not deviate very far from its desired depth. Plots with the same pair of 
coefficients (ADE and ADR), with the same dynamic length of cable, and 
with the same type of inital condition perturbation, but with different 
values of spring constant, seem to be very similar. The responses with 
the larger spring constant (SC=1500) seem to settle more quickly, even 
though they are oscillating through a wider range of depths. All the 
responses to the fin angle perturbations are remarkably similar. This 
might be because the fin is quickly driven back to the desired stability 
position; and since all the fin angle perturbations are the same pair of 
values, this similarity in responses is not Surprising. The perturbing of 
the initial conditions for the new tension force realization was the point 


at which this study ended. 


ADE=1.2 
ADR=0.07 
APR=3. 
SC=500 
All -2.,1—+0,, 9 
All 2 aU 
FAI -17.5-+15.0 
SC=1000 
Z1I te ea 
All ale, 5 + Que 
FAI ie ie) Saal) 
SC=1500 
Z1] -4.8-43.1 
All —2Oe FLU. U 
FAI =F .0-Fl9.0 


ZO) = 5 ll 


ADE=1 .2 
ADR=0.3 
APR=3 , 


me oT Ue 


+3 ,.0-+7.5 


aiveso— +150 


— 2D: Aaanl 


Opi: Ora, ..0 


-17 .o- +150 


-2 7743.3 


A Ua aS 


-17.5-4+15.0 


Table V 


ZQO=100., 
ii). 2 ADE=1..2 
ADR=0.07 ADR=3 .0 
APR=3. APR=3. 
-1.7-+1. -1.5-40.9 
+2 .0-+7. +2 ..5->4+6.5 
-17.5~415.0 -17.5-4+15.0 
-4,3342.5 -2.8742,2 
-0.5-+8.0 0.0-+8.0 
-17 .5-415.0 -17.5~415.0 
-4,8-+3.0 -2.7-42.9 
-1.5--+9.0 -1,0-+4+8.5 
-17.5-+15.0 -17.5-4+15.0 


STABILITY RANGE OF INITIAL CONDITIONS 


V. CONCLUSIONS AND RECOMMENDATIONS 

The study presented here is only a beginning to this problem. There 
are many facets and additions to this problem that remain to be Studied. 
First of all, the dynamic characteristics and wave-following capabilities 
of this model were never investigated. Other additions could be made to 
the model to take into account the heaving forces generated by the waves 
and the free-surface effect which will definitely affect the motion and 
capabilities of the buoy. There remains much to be studied in the area 
of the dynamic cable and the coupling of the tension force to the buoy 
model. Other tension force realizations might be tried with several spring 
constants in a single cable length; or the cable might be approximated by 
a series of straight lengths, each with its own spring constant. The 
possibilities of adding a damping coefficient into the tension force real- 
ization need to be investigated; the possibilities are almost endles. 
Other methods of simulating the buoy need to be tried, especially that 
of using a hybrid computer (an analog computer integrated with a digital 
computer); hybrid computers are more suited to this type of real-time 
Simulation work. A dynamic model that allowed for all the possible 
dynamic forces and all the possible variations in buoys and cables, 


would be invaluable in the design and engineering of towed buoys. 


ay 


APPENDIX A 
RESPONSES OF THE BUOY TO INITIAL CONDITION PERTURBATIONS 
The following 24 plots are the responses of the buoy to the maximum 
initial condition perturbations listed in Table V for SC=500 and SC=1500. 
The plots are of depth versus time. The x scale for all the plots is 


5 sec/inch; and the y scales are in ft/inch. 
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APPENDIX B 


SYMBOLS AND ABBREVIATIONS USED IN THE COMPUTER PROGRAMS 


Z1 


is 


Py 


Pal 


Xl 


XZ 


PH 


Te 


ala 


PG 


PE 


SFE 


CFE 


ior 


Cr 


VAN 


V2 


vertical position in ft. (initial condition ZI) 
vertical velocity in ft/sec. (initial condition Z11) 
vertical acceleration in ft/sec* 

Summation of vertical forces 


Summation of the vertical forces less the vertical tension 
Component 


horizontal velocity in ft/sec. (initial condition X1I) 
horizontal acceleration in ft/sec* 
Summation of the horizontal forces 


summation of the horizontal forces less the horizontal 
tension force 


the absolute value of Fl 
the absolute value of F2 


the angle with which the tension force acts on the body 
in radians 


the sine of FE 

the cosine of FE 

the vertical component of the tension force 
the horizontal component of the tension force 
horizontal velocity in knots 


the angle of the velocity vector relative to the horizontal 
in radians 


the magnitude of the velocity vector squared 


Al 


Yel 


WZ 


YM 


FA 


FB 


Ae 


AB, 


BR 


CL 


CD 


CM 


CV 


OV 


AD 


the body angle of attack relative to the horizontal in radians 
(initial condition AlI) 


the body pitch rate about the Y axis in radians/sec 
(initial condition Y1I) 


angular acceleration of the body about the Y axis in 
radians /second* 


summation of the moments on the body about the Y axis 


the body angle of attack relative to the velocity vector 
in radians 


AR in degrees 


the fin angle of attack relative to the horizontal in 
radians (initial condition FAI) 


used in the fin servo loop 
used in the fin servo loop 


summation of the inputs and the feedback in the fin 
servo loop 


the fin angle of attack relative to the velocity vector in 
radians 


BR in degrees 

coefficient of lift of the body 
coefficient of drag of the body 
coefficient of moment of the body 
the temsiionsience 

the cosine of VAN 

the sine of VAN 

the absolute value of z 


the antenna drag 
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BL 


BD 


YCM 


YADS 


YADC 


ZD 


DEE 


DE 


DRR 


DR 


YPR 


ADE 


ADR 


APR 


the body lift 
the body drag 
the moment on the body about the Y axis 


the moment on the body due to the vertical component of 
the antenna drag 


the moment on the body due to the horizontal component 
of the antenna drag 


in the first program this is used to measure the relative 
depth error from 4.0 ft; in the second program it is used to 
compute the antenna drag 

the depth error in psi/ft 

DEE times the coefficient ADE 

the depth rate in psi-sec/ft 

DRR times the coefficient ADR 

the body pitch rate (Yl) times the coefficient APR 

the depth error coefficient 


the depth rate coefficient 


the pitch rate coefficient 


THE FOLLOWING ARE USED SOLELY IN THE SECOND COMPUTER PROGRAM 


PSS 


ZO 


SC 


OFS 


LO 


IN THE REALIZATION OF THE NEW TENSION FORCE 
the steady-state value of FE 


the vertical approximation to the dynamic length of the 
cable; see Fig. 17 


the spring-constant value of the cable in pounds/ft 
the sine of FSS 


the approximation to the dynamic length of the cable; 
see Fig. 17 


TE 


XO 


ZN 


XN 


LN 


FE 


the tangent of FSS 


the horizontal approximation to the dynamic length of the 
cable; see Fig. 17 


the new computed value of ZO 


the relative movement out of steady state in the X 
direction 


the change in position out of the steady state in the 
X direction 


the new computed value of XO 


the new computed length of the dynamic length of the 
cable 


the change in length of the cable 


as stated previously, but computed differently 
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AN ADDITIONAL SET OF PARAMETERS TO BE STUDIED 
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GRAM CONTROL CARDS 
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CON 
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END 
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STOP 


NOT INCLUDED IN ACTUAL PROGRAM 
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